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Abstract 

In order for finding a good individual for a given fitness function in the context of evolutionary computing, we 
introduce a novel semiclassical quantum genetic algorithm. It has both of quantum crossover and quantum mutation 
procedures unlike conventional quantum genetic algorithms. A complexity analysis shows a certain improvement 
over its classical counterpart. 

1 Introduction 

■ Continuous development has been performed on genetic algorithms [iHSj. Along with the development of quan tum 
\ computing [1,0, quantum-inspired classical algorithms for evolutionary computing have been developed (6l4l3l| (see 
^ • also a review [IJ] and references therein ) . In addition, classical genetic algorithms to evolve quantum circuits have also 



been studied by several authors |15l - l29l | (see also review articles (30|,|31|). These algorithms are, however, designed to 
work on classical computers. Quantum genetic algorithms (QGAs), in its literal meaning, nonetheless, have gathered 
comparably little attention and a few works |32h36| | have been performed so far. Evolutionary computing on quantum 
• architectures will achieve more attention if there is a scenario to establish significant improvement over classical 
[ counterparts. Indeed, Malossini et al. H, 2^ claimed that computational complexity of an evolutionary step between 



■ generations is exponentially fast in their quantum algorithm in comparison to a classical one. It seems, however, that 
C^l ', they overlooked the complexity of index-to-string conversion circuits or, otherwise, overlooked the circuit complexity 
' of a variant of the inversion-about-average operation. Their claim is thus misleading as we will discuss in Sec. [D 



] Recently, Johannsen et al. 37[ applied quantum search algorithms to several optimization problems in a certain 
OA ■ context of evolutionary computing. Nevertheless, it is obscure how much cost is spent for the internal quantum circuit 
of a variant of the inversion-about-average operation used for amplitude amplification (namely, the operation denoted 



K*" ■ as ^Sq^~^ in the convention of Ref. [38|) in their approach. 

k>( , Let us briefly summarize the conventional approach of QGAs and its problem. Here, we omit discussions on 
' incomplete or physically unfeasible works on QGAs, which were summarized by Sofge in Ref. (39| . The aim of genetic 
Ci \ algorithms is, in most cases, to find an individual (typically, an input string) with a very high fitness value for a 
given problem. It starts with 0(1^1"") initial individuals where \A\ is the size of the alphabet A and w is the length 
of a schema expected to be a building block for a given problem. Let us restrict the problem by representing each 
individual as a chromosome encoded as a binary string with the length n. We need ^ N — 2^ initial chromosomes 
for the worst case. With a quantum register, one may use n qubits to make a superposition of N chromosomes. We 
regard the probability to find a chromosome on the computational basis as its (normalized) population. 

First we briefly overview the selection strategies in conventional QGAs. In short, a selection is an operation to 
enhance populations of individuals with high fitness values and to decrease those of low fitness values. In quantum 
computing, it is natural to utilize variants of Grover's algorithm [40| for this purpose. In fact, a variant of the Grover 
search for a priori unknown number of solutions (the Grover-BBHT search) [4l| and that for finding the maximum 



42, |43| are the essential parts of the QGAs developed in Refs. |33l436| . The query complexity 0{\jN/r) for a variant 



of the Grover search dominates the total complexity in the QGA of Udrescu et al. [33|, |35| . For the QGA of Malossini 



et al. 3J, |36j, there is a different factor to consume time. As mentioned, this will be discussed in detail in Sec. [2] 

Second we overview the strategies for crossovers and mutations in conventional QGAs. In a short explanation 
of these terms, a crossover is to exchange substrings of two chromosomes and a mutation is to flip certain bits of a 
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chromosome. These are effective operations to enlarge the search space. So far, quantum crossover operations have 



not been developed. Malossini et al. [3J, [Sg used classical crossover and mutation operations; Udrescu et al 
did not use crossover and mutation operations. Johannsen et al. [ST*! introduced quantum mutation operations in 
an application of quantum amplitude amplification [38.] to some optimization problems; nevertheless, crossovers were 
not used. It is, in fact, in general difficult to manipulate populations of quantum states for handling crossovers if we 
imitate the classical way in a straightforward manner (See also Ref. js^l which explained the difficulty in a slightly 
different manner). Picking- up two particular individuals and making a crossover costs 0(n) quantum gates if we use 
0{n) ancillary qubits. Since there are many possible pairs for a crossover at each generation, there is no speedup 
over a classical crossover. More specifically speaking, this approach needs to look up classical data of chromosomes 
to specify a pair of chromosomes. Thus there cannot be any speedup. In addition, this requires an exponentially 
large classical memory in comparison to the size of a quantum register. We will show a different approach to handle 
crossovers in our algorithm in Sec. |31 

As another direction for developing quantum evolutionary cornputing, which we do not pursue in this contribution, 
one may use so-called quantum fixed-point search algorithms 3, • They increase the population of a superposition 
of target chromosomes by iterative applications of unitary amplitude-amplification operations. Unlike the original 
Grover search, there is no drawback in an excessive iteration; the probability of finding a target grows monotonically 
as 1 — e^'+i where e is the probability of finding a non-target chromosome and t is the number of queries. The query 
complexity of this approach is, however, as large as that of the exhaustive search: When the number r of targets is 
very small in comparison to iV, we have e^*"*"^ = (1 — r/NY*'^^ ~ 1 — (2i -I- l)r/N, which implies that 0{N/r) queries 
are required to achieve a sufficiently small error probability, say 1/2. Thus for r ^ N, Grover's algorithm should be 
chosen instead of fixed-point quantum search algorithms. In addition, the asymptotic optimality of e^'+i as an error 
reduction speed was proved [46| for fixed-point quantum search methods. It is also known that the query complexity 
of the Grover search is optimal 41 1 in general as a unitary process for unsorted search. Therefore, it is unlikeljo that 
a fixed-point quantum search is effectively used in an evolutionary computing instead of the standard Grover search 
and its variants. 

In this contribution, we propose a semiclassical quantum genetic algorithm that involves a quantum crossover 
and a quantum mutation. It uses the quantum search for finding the maximum 42. |43| for the selection procedure. 
Although this selection strategy looks similar to those of conventional QGAs 33"36|, it is different in the point that 
we use tailored inversion operations for the quantum search. Our algorithm achieves a quadratic speedup over its 
classical counterpart as an algorithm involving all possible crossovers for a chosen crossing site. 



The paper is organized in the following way. Section [2] discusses on the algorithm of Malossini et al. [3J, |36| . 
Section [3] describes our algorithm. The procedures of the crossover, the mutation, and the selection are developed 
in the standard quantum circuit model. Computational complexities of the algorithm are evaluated in Sec. SI The 
obtained results are discussed in Sec. [S] and summarized in Section [S) 



2 A conventional quantum genetic algorithm 

Here we discuss the conventional quantum genetic algorithm proposed by Malossini et al. js^, [3^ • There is another 
conventional algorithm proposed by Udrescu et al. 33, 3^, which is not introduced here. The main difference is that 
Udrescu et al. did not employ crossovers and mutations while Malossini et al. employed them as classical operations. 

The problem we consider is given as follows. 
Problem: Suppose there are N individuals with indices 0, . . . , iV — 1. There is a fitness function / : {0, . . . , — 1} ^ 
[0, 1]. Find one of the individuals with sufficiently large fitness values. 

The algorithm of Malossini et al. is described as follows. 



Algorithm 1 (Malossini et al. [3J, |36|) 

Main Procedure: 
Start from a generation consisting oi N < N individuals. 
Repeat: 

1. Use the Quantum Selection Procedure to get one index. Do the same procedure to get another. 

2. Classically make a crossover for the two individuals corresponding to the indices obtained in 1. Clas- 
sically make a proper mutation. 

3. Classically replace randomly chosen two individuals with the two offsprings obtained in 2. 



^ There is another drawback in a quantum fixed-point search: it requires a nonstandard oracle operation represented by the unitary 
gate Uo = l — 2\t){t\ with \t) a superposition of target states. A standard oracle is represented by the unitary gate Uo = 1 — 2 \Tm){Tm\ 
with \Tm) the target states. There is no known way to construct Uo from Uo efficiently when we have multiple targets that are unknown 
in advance. This difficulty exists not merely in the fixed-point search but also in a variant of Grover search using Uo [47| . 
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Quantum Selection Procedure: 

(i) Choose randomly an index y among N possible ones in the present generation and compute its fitness 
/(y). Set the threshold Fy i — /(y). 

(ii) Perform t times: 

(a) Initialize the quantum registers to — ^ l^^)!?/)! where x's are the indices of the present gener- 

V N 

ation. 

(b) Apply the Grover-BBHT search |4l|, where we internally use an oracle that inverts the signs 
of the marked states, namely states \x) such that f{x) > Fy. This enhances the amplitudes of 
marked states after its Grover iteration. 

(c) Measure the left register in the computational basis and get the new index y'. If f{y') > Fy, then 
yi~y' and Fy ^ f{y'). 

(iii) Return the index y. 

They analyzed their algorithm and claimed that the selection procedure internally requires a very small number of 



oracle calls in contrast to 0{N log N) oracle calls required by a classical selection procedure [3J, |36|. In short, their 
claim was that the number of the Grover iterations in step (b) was 0(1). 

Their complexity analysis was, however, misleading. A more accurate description of the complexity is given in the 
following way. There are two possible cases we should consider: (A) TV is approximately equal to N; (B) N ^ N holds. 

Let us consider the case (A). We can preset the index x to represent the a;th individual in a good approximation 
thanks to N c:i N. Then the given fitness function can be used as it is, without index conversion. This indicates that 
we use nothing but the standard Grover-BBHT search for N indices. As TV ~ TV, the quantum selection procedure is 
called only once or a very small number of times. Suppose there are r nearly-optimal individuals that are acceptable 
as solutions for a given problem. Inside the quantum selection procedure, as the selection goes on in (ii), the oracle 
used in step (b) has less marked states. The number of marked states converges to r. Therefore, the query complexity 



of a single step in a very later selection stage is 0{y N/r) = N /r). This is the accurate description of the query 
complexity for the present case and it is same as that Udrescu et al. gave for their algorithm 3^ 3^. On average. 



each time of the repetition in (ii) virtually extinguishes a half of the individuals by increasing the threshold. Thus, by 
considering the sum of the geometric series N/{2^r) {i = 0, 1, 2, . . .), the average total query complexity is found to 
be 0(a/]v7^). 

Let us now consider the case (B). In this case, the quantum selection procedure is called several (or more) times. 
The important fact is that the chromosomes of individuals in the generation are not consecutive binary strings. The 
generation contains N binary strings out of N possible ones. There are two options to handle them: (B-1) consecutive 
integers are assigned as pointers to nonconsecutive chromosomes; (B-2) nonconsecutive chromosomes are used as they 
are. 

(B-1): In this case, the index x is a pointer to the chromosome Kx for which the fitness function works appropriately. (In 

short, this is a workaround for the fact that the original Grover-BBHT search should start from a uniform superposition 

of consecutive indices while the generation consists of nonconsecutive chromosomes.) Therefore, for the superposition 

\iI;q) — —j= ^ |a;)|j/), the fitness function cannot be applied directly. We need a conversion circuit to interpret x as 
V N 

Kx- Then the fitness function is applied and it returns f{nx)- Each conversion should involve 0(log A^) CCNOT gates. 
Thus, the total circuit depth for the index conversions is 0(A^log A^). As a result, a single oracle call accompanies the 
additional time complexity 0(A^log A^). Indeed, the query complexity is small if N is much smaller than N. However, 
the cost of index conversions hinders speedup over a classical selection procedure. 

(B-2): In this case, the state \x) in the superposition |-0o) — —j^ keeps the actual chromosome x as an index. 

V N 

One should use the Grover-search routine somehow without index conversions by choosing |i/;o) as its initial state. 
Indeed, this is possible [s^l in the case where one can provide the operation L = I — 2|-0o)(V'o| as the inversion-about- 
average operation instead of the standard one for the original Grover search. Nevertheless, L should be generated by 
sandwiching / — 2|0)(0| with J7(a) and Uj^^^ in the present context, where the unitary operation [/(a) corresponds to the 

initiahzation step (a). Thus the circuit complexity of L is 0(A^log A^) (see, e.g., Ref. [il]). Such an expensive circuit 
is used subsequent to every query. Obviously, the time complexity in this case is as large as the one in case (B-1). 

One may also think of replacing the Grover-search routine with the generalized Grover search for nonconsecutive 
integer sets (49| . Nevertheless, it is required to find an appropriate unitary transformation replacing the standard 
Hadamard transformation H ® ■ ■ ■ ® H . There is no known way to find it efficiently when a random integer set is given 
and target integers are unknown in advance. 

There is, in fact, a way to reduce the circuit complexity of L introduced in (B-2) if the algorithm is modified so 
that it uses an efficient pseudo-randomizer (or, pseudo-scrambler) circuit instead of directly using a random number 
generator. Our algorithm introduced in the next section takes this approach. 
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In addition to the above discussions, we should mention that Algorithm[l]has another problem: it takes 0{N log N) 
space to keep the generation in a classical memory. This is usually quite larger than the 0(log A^) space that is enough 
for quantum search (neglecting the space internally used by an oracle circuit). Our algorithm is designed not to face 
this problem. 

3 Algorithm with quantum crossover and quantum mutation 

Conventional quantum genetic algorithms have not achieved improvement in the total complexity over the standard 
quantum search, and quantum crossover operations have not been developed so far, as we have seen in the previous 
sections. Here, we show an algorithm using a quantum crossover and a quantum mutation. We also use a variant of 
the Grover search for the selection procedure, whose query complexity is the main factor of the total time complexity. 
We have still a certain speedup over a classical counterpart as a genetic algorithm using a crossover. It is novel in 
the sense that a simultaneous crossover using a superposition is achieved. The speedup partly relies on an efficient 
internal structure of a pseudo-randomizer circuit, which will be explained in Appendix. 

3.1 Algorithm flow 

We introduce our algorithm as follows. Here, the quantum register is accessible from subroutines as a kind of global 
variables. 

Algorithm 2 

Consider a threshold /th to the fitness values, which is considered to be sufficiently large. 
t i — 0. 

REPEAT 1.-9.: 

1. Construct a pseudo-randomizer R that maps a c-bit string to a pseudo-random n-bit string. It 
should be implemented as a circuit whose input and output are integer couples (a,0) and (a, i?(a)), 
respectively, where a is a c-bit integer. The circuit should consist of poly(cn) elementary reversible 
logic gates. We require R to regard Oq • • • Oc^i as an exception and map it to Oq • • • 0„_i. We also 
require N = 2'^ ^ N — 2"-. Once constructed, R is fixed until next t. An explicit example to construct 
R is given in Appendix. 

2. IF t = THEN generate a random c-bit string 7 and set z i — R{'j) (otherwise, z is the best 
chromosome found in the {t — l)th trial) ENDIF 

3. Generate a random c-bit string 7' and set u i — Ril')- 

4. CALL init_reg(i?, z) defined in Sec. 13.21 twice to make two identical quantum states that are both 
\(p) = — t^X^kpx ka:)a|2^) whcrc is the address pointing to x £ X; X = {i?(g)}„ U {z} with 

g = Oo • • • 0c-2lc-i, • • • , lo • • • Ic-i; subscript "a" stands for the address portion. We write the entire 
unitary operation of this procedure as C/init- That is, \(p)^^ — ?/init(|0)a|0))'^^. 

5. CALL quELntuin_crossover(/) defined in Sec. 13.31 for the current quantum register, which is a 1-point 
crossover with the crossing site placed between the {I — l)th and the Ith bits for each chromosome. 
This procedure is an identity map as a quantum operation acting on \lp)'^'^. 

6. Apply the quantum mutation (Sec. 13. 4p to the current quantum register. We write the entire unitary 
operation of this procedure as t/mut- Apply the same mutation classically to u. 

7. CALL quantmn_selectioii(J7init, CAnut, u) defined in Sec. 13.51 for the current quantum register and 
obtain the output chromosome z. 

8. IF /(z) > /th THEN RETURN z and EXIT ENDIF 

9. t i — t + 1. Refresh the quantum register. 

3.2 Preparing the initial state of a quantum register 

In this section, we define the procedure init_reg(i?, z). 

This procedure is intended to prepare a superposition corresponding to the generation given as a set X of n-bit 
integers x with jJA" = N — 2'^. The procedure starts with the quantum state |0)f"^|0)®". The desired superposition 
is \ip) = —j= ^ px \^x}a\x) with Qx the address pointing x. We opt to use consecutive addresses 0, . . . , 2^^ — 1. The 

V -'V 

order of time complexity of the procedure equals to the internal circuit complexity poly (en) of the pseudo-randomizer 
R. An explicit example to construct i? as a quantum circuit is given in Appendix. 
The procedure is defined as follows. 
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PROCEDURE init_reg(X): 

(i) make a superposition ^ Ejlo^ b')a|0 • • • 0) by applying H^" ® to |0)f =|0)®". 

(ii) Apply the pseudo-randomizer R implemented as a unitary operation mapping |j)a|0) to with 
Xj an n-bit pseudo-random number for j = 1, . . . , 2*^ — 1. 

By assumption, |0)a|0) is mapped to |0)a|0). 

(iii) Apply a Oq • • • Oc-i-controUed X^" (g) • • • (g) to map |0)a|0) to |0)a|^;), where Zk is the kth bit of 
z (fc = 0, . . . , n — 1). 

(iv) RETURN the current state, namely \ip). 



3.3 Crossover operation 

Here, we construct a quantum crossover procedure quaiitum_crossover(Z). It is a 1-point crossover acting on the 
chromosomes simultaneously. 

Recall that the original generation is given as a set X of n-bit integers x with {jX = N. Each x has its left side 
x'''" and its right side x"^^*^ separated by the crossing site. The crossing site is placed between the (/ — l)th and 
the Ith qubits as specified by the parameter. Hence, the bit length of a;'*''* is I. With the procedure, we generate a 
superposition of all the children that are combinations of x^''**'s and x"^^^'s together with their parents with the same 
weight as children. 

The state of the quantum register in the beginning of this procedure is 

This state has the components |a;^'^**)|a;'''s'^*)|a;''''^*)|a;''^'^'^*) besides the addresses. We relabel the qubits so that the 

middle portion 

l^right^l^/left^ is put aside from our minds. Let us conceal them by denoting as l^xa;')- addition. We 
denote the main portion jaj'^^'x'"^'* ) with the subscript "main". The state |<^)®^ with the new qubit labels is written 
as 

— V V If? O Uleft„/"ght, I , 

We have at most (ftX)^ distinct chromosomes in this state. In this way, all possible crossovers are performed at once by 
the relabelling. The resultant state is a superposition of all of the children together with their parents. (The parents 
are involved because the values of x and x' may coincide.) This is desirable as a crossover because sometimes some 
parents have higher fitness values than any child. 

The procedure described above is formally written in the following way. 

PROCEDURE quantum_crossover(Z): 

(i) We have the quantum register in the state The original labels of its qubits are 0, . . . , 2c + 

2n — 1. We relabel them as 0, . . . , c — 1, 2c, . . . , 2c -|- Z — 1, 2c + n, . . . , 2c -|- 2n — Z — 1, c, . . . , 2c — 1, 




(ii) RETURN 

As is obvious, this procedure is an identity map as a quantum operation. Once the crossover process is completed, 
one may use a mutation as an option. This is going to be explained in the next section. 



3.4 Mutation 

In classical genetic algorithms, randomly selected chromosomes are affected by a mutation, such as bit flips acting on 
randomly-chosen places. Here, we consider the mutation described below. Although it is written as a classical routine, 
it can be trivially interpreted to a quantum circuit. 

PROCEDURE tinp_inut: 

(i) Let us randomly generate the first template like ***0*1***1**0****1*, which specifics a schema to 
mutate. Using this template, we pick up chromosomes with specified bits like 0,1,1,0,1 in the specified 
places. 
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(ii) We also use the second template like *x*************X*** in which X's can be placed only on the 
places where *'s (namely, "don't care" symbols) are placed in the first template. Using this template, 
we apply the bit flip X to the specified places of the chromosomes picked up in (i). 

(iii) RETURN 

Note: Technically, we often wish to avoid a mutation for the best chromosome z found so far by the present 
time step. This is realized by choosing the first template so that this does not happen. 

As a quantum circuit, tmp_mut is realized by a multiple-bit controlled multiple-bit NOT gate. In the present example, 
the gate is "0-controlled 1-controlled 1-controlled 0-controlled 1-controlled NOT NOT" with control bits 
specified by the first template and the target bits specified by the second template. The gate acts on the portion 
l^icft^/nght^^^^^. the addresses and ax' are untouched. Therefore, the resultant state can be written as 

, ^, — right 

where (x x' )'s are the chromosomes after the mutation process. 

The next step is to apply a natural selection to the chromosomes living in the superposition 



3.5 Selection 

In this section, we introduce our selection procedure. It is intended to find a chromosome having the maximum 

fitness among (a;'°^'a;' ^ )'s. It utilizes the quantum search for finding the maximum 42, 4^. As we have mentioned, 
conventional QGAs |33l436l ] have similar selection strategies. The difference from them is that we use tailored inversion 
operations for the quantum search. 

This procedure is called with three arguments: J7init, C^mut, and u. We have the state ([T|) at the beginning of this 
procedure. The procedure is formally defined as below. Here, /ctcrm = i] x [(45/2)A^+ (28/5)(log2 A^)^] with integer 
constant t] > 1. 

PROCEDURE quantum_selection(t/init, C^mut, u): 
FOR ki — TO fctorm - 1: 



(i) Set 



Ui=I®I®I-2 V I®\y){y\®I 



f{y)>fiu) 



where the left and the right Fs act on the address states {a^ax') and the states \*xx'), respectively. 
This is the oracle function that can be implemented as follows. First, we attach ancillary qubits as 

blocks (I) and (II) in the state |0 ••• 0)(i) |-)(n) with]-) = (|0) - |l))/\/2. For each xx', /(J^<=f*a/"^^') is 
set as a value of the block (I) by using / implemented as a quantum circuit. Let us write this operation 

as Uf. Then we compare the block (I) with f{u); we fiip the qubit (II) if /(a;'°^*x' * ) > f{u). We 
apply Uf^ (this disentangles the ancillary qubits from the main register) and remove the ancillary 
qubits. 

(ii) Set 

U2 = /(g)/®/- 2|<^)(^|. 
This operation is composed in the following way. 

U2 = C/mutC/i„it[/ - 2(|0)«^|0)^"(0|«^(0|^")®2][/t^.^C/t^^^, 

where we also relabel qubits according to quantum_crossover(^). 

(iii) Apply the Grover-BBHT search [4ll| for which we use Ui and U2 instead of the standard operations, 
namely, the inversion operation for targets and the inversion-about-average operation, respectively. 
As for the starting state of the search, we use \(p). 

(iv) Measure the main register and obtain the chromosome u' . 
IF f{u') > f{u) THEN u ■ 

END FOR 
RETURN u. 
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Let us explain on the number of iterations, ktcm- The above procedure is same as the weh- known quantum search 



algorithm for finding the maximum 4^, |43[ except for the definitions of the inversion operation for targets and the 



inversion-about-average operation. In other words, we perform the quantum search for finding the maximum in the 

subspace spaii{\aa:ax')a\x^'^^^x' ^ )main|*a;a;')}- As provcd by Diirr and H0yer [i^], the probability for the output to 
be the maximum is at least 1/2 if the number of iterations is [(45/2)-\/M -|- (7/5)(log2 M)^] with M the number 
of indices. In the present context, we have M — since there are iV^ distinct addresses. Therefore, after fctcmi 

, r. — right 

iterations, we find the output chromosome having the maximum fitness among [x x' )'s, with the probability at 
least 1 — (1/2)''. By setting -q sufficiently large, say, around 16 to 24, we have the desired output with almost certainty. 



4 Computational cost 

We evaluate the computational cost of each process of Algorithm [2] in Sec. 14.11 The total computational cost to handle 
a generation is derived in Sec. 14.21 and compared with that of classical counterpart in Sec. 14.31 

4.1 Costs of each procedure 

For a single call of each procedure, the costs are evaluated as follows. 



Cost of the initialization The procedure to prepare the initial state, described in Sec. 13. 2[ uses poly(cn) = 
poly(log A^log A^) elementary quantum gates. 

Cost of the crossover The crossover described in Sec. 13.31 makes use of two identical (c + n)-qubit states. This 
procedure does not use any quantum operation but relabels qubits. This takes 0(log A^) time. 

Cost of the mutation The mutation described in Sec. 13.41 involves a single gate that looks like, say, 
* - NOT - Co - *-Ci - Ci-*-* - Co-*-* - NOT - * - Ci - * with multiple control bits (0- 

control Co's and 1-control Ci's) and multiple NOT gates placed according to the corresponding templates, where 
symbol * stands for an untouched qubit. Such a gate can be realized by O(logiV) elementary quantum gates with 
0(log A^) ancillary qubits. 

Cost of the natural selection The natural selection described in Sec. 13.51 consumes 0{N) queries to find the best 
chromosome among the set of at most N'^ chromosomes. Each query accompanies the operation U2 that invokes C/init 
and C/mut, and also their inverse operations. It is easy to find that the internal cost of U2 is poly(cn) according to 
the costs for C/init and t/mut- In addition, the cost to prepare the starting state is poly (cn). We may also mention 
that, the internal cost of the fitness function is a certain small factor, typically poly(n), as conventionally assumed Q. 
Therefore, the circuit complexit}]! of the procedure is 0(7Vpoly(log7Vlog A^)). As for space, we use poly(log A^) qubits 
in total, considering a typical fitness function consuming poly(n) space. When the fitness function is designed to use 
0{n) space, 0(log A^) qubits are enough, although we do not assume this case for evaluating the space complexity. 



4.2 Total cost 

Comparing the costs of the four procedures, the dominant cost is the circuit complexity for the selection procedure. 
Therefore, we find that our algorithm uses 

0(7Vpoly(log7VlogAr)) (2) 

elementary quantum gates for each t. This is the time complexity of our algorithm for handling each generation. 
As for the space complexity, we spend poly (log A^) qubits and O(logA^) classical bits as is easily evaluated from the 
description of Algorithm [21 



4.3 Comparison with a classical counterpart 

The classical counterpart of our algorithm is the following algorithm. As is obvious from its structure, it should have 
the same output and the same number of iterations as our algorithm with almost certainty as long as the same random 
seed is used for step 1 to construct R for each value of t. 

Algorithm 3 



■^We assume that the circuit depth is on the order of the circuit size, namely, the number of elementary quantum gates. 



7 



Consider a threshold fth to fitness values, which is considered to be sufficiently large. 
ti — 0. 

REPEAT 1.-7.: 

1. Construct the pseudo-randomizer R found in Algorithm [2] as a classical function. Fix R until next t. 
Using R, we generate a generation X = {R{q)} with q ~ Oq ■ ■ ■ 0c-2lc-i, • • ■ , lo • • • Ic-i- We require 
N = 2" = 2\ 

(Thus, we have A — 1 ?i-bit-lcngth chromosomes in X presently.) 

2. IF t 7^ THEN put z into X {z is the best chromosome found in the (t — l)th trial) ELSE generate 
a random c-bit string 7 and put R{j) into X ENDIF 

(Now we have A chromosomes in A.) 

3. Split all cc € A between the {l~ l)th bit and the Ith bit. This makes Z-bit strings x^'^^^^s and (n— Z)-bit 
strings a;'''s^''s. Generate the set A consisting of all x's and all of their children that are all the 
combinations of a;'''^*'s and a:'^'sht'g_ 

4. Apply the template-based mutation tmp_mut, introduced in Sec. 13.41 as a classical procedure to all 
chromosomes in A. 

5. Find the chromosome z having the best fitness value among those in A. 

6. IF f{z) > fth THEN RETURN z and EXIT ENDIF 

7. ti — t + 1. 

The computational cost of each procedure of Algorithm |3] is as follows. 

• Generating A in the steps 1.-2. takes 0(Apoly(cn)) basic operations. 

• We need to use 0{nN^) space and 0{nN^) basic operations to perform all the crossovers among A parents. 

• Mutations acting on the individuals of A takes 0{nN^) basic operations. 

• The selection to find the best individual from A takes 0(A^poly(n)) basic operations considering the cost 
poly(n) of calculating a fitness value. 

As n = log2 A and c = log2 A, the time and space complexities are 0(A^poly(log A)) and 0{N^ log A), respectively, 
for each t, i.e., for handling each generation. 

In contrast, as we have seen in Sec. 14.21 the time and space complexities of our semiclassical quantum genetic 
algorithm (Algorithmic]) are 0(Apoly(log A log A)) and poly(logA), respectively, for each t. Therefore, neglecting 
the difference between poly(log Alog A) and poly (log A), we have achieved a quadratic speedup over its classical 
counterpart together with an exponential reduction in space. 

The classical counterpart has been constructed by keeping the one-by-one correspondence with the quantum al- 
gorithm. Thus there is a possibility that a better classical algorithm with the same behavior exists. This is in fact 
the case for Algorithm [3] For a fairer comparison, now we reform Algorithm |3] and reduce its space complexity. The 
following algorithm has the same output and the same number of repetitions as Algorithm |3] while its space complex- 
ity is exponentially reduced. We use the same pseudo-randomizer construction and the same mutation procedure as 
before. Of course, we set N — 2'^ and A = 2" for integers c and n satisfying 1 < c < n. 

Algorithm 4 

Consider a threshold fth to fitness values, which is considered to be sufRciently large. 
t i — 0. 

REPEAT 1. 8.: 

1. Construct the pseudo-randomizer R : {0, 1}'^ — {0, 1}" as a classical operation. 

2. Construct the mutation process as a map M as a classical operation. 

3. IF i = THEN for a random 76 {1, . . . , A — 1}, z < — ^(7) (otherwise, z is the best chromosome 
found in the (t - l)th step) ENDIF 

4. For a random 7' e {1, . . . , A - 1}, j < — R{l')- 

5. FOR a< — TO A - 1: 

IF a ^ THEN xi — z ELSE x i — R{a) ENDIF 
FOR b< TO A - 1: 

IF & = THEN yi — z ELSE y i — R{b) ENDIF 

Crossover x and y and obtain children v and w. 

Find the best chromosome g among the chromosomes M{x), M{y), M{v), and M{w). 
IF fig) > f{j) THEN J ^ g ENDIF 
END FOR 
END FOR 
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6. z i — j. 

7. IF /(z) > /th THEN RETURN z and EXIT ENDIF 

8. ti — t + 1. 

In this algorithm, R is called 0{N'^) times for each value of t. We know that R internally takes poly(cn) time. Therefore, 
this algorithm spends 0(A^^poly(log N log N)) time for handling each generation. As for the space complexity, it spends 
only poly(log N) space, which is clear from the algorithm structure. 

In comparison to this enhanced classical algorithm, our semiclassical quantum algorithm still has a quadratically 
small time complexity as shown in Eq. ([2]). 



5 Discussion 



How to perform crossovers in a quantum manner was a pending problem in conventional quantum genetic algorithms 



33H36|. In fact, a selective crossover for specific two chromosomes is expensive when they are component states of 
a superposition. Even if we attach address states pointing to the component states, we need to look up the classical 
data of the chromosomes to construct a quantum circuit realizing the unitary operation for this purpose, or more 
specifically, to place address-controlled bit-flips appropriately. In our algorithm (Algorithmic]), we have avoided to 
mimic a classical way and chosen a different approach. We use two identical copies of a superposition corresponding 
to a generation and utilize relabelling of qubits so as to handle all possible combinations of substrings simultaneously. 
Obviously, the classical counterpart of our algorithm is the one that seeks for the best chromosome (after a mutation) 
among all possible crossovers for a chosen crossing site for each generation. Comparing our algorithm with the classical 
counterpart, we concluded that we have achieved a considerable reduction in the computational cost. 

One may, however, claim that usually at most several crossovers are performed for a single generation in a classical 
genetic algorithm. Indeed, our algorithm is not aimed to be a quantum counterpart of a common classical genetic 
algorithm. As we discussed in Sec. [21 a straight-forward conversion of a common classical algorithm into quantum 
one by simply incorporating a quantum search into the selection procedure has a problem: we need to either interpret 
nonconsecutive integers to consecutive ones or use an expensive construction for the inversion-about-average operation, 
in order to perform the Grover-BBHT search, which causes a significant loss of performance. This problem should be 
resolved so as to find an appropriate quantum counterpart for the common case. Seemingly, the following work-around 
looks fine: (i) Use the pseudo-randomizer R used in Algorithm [2] instead of a random number generator to generate 
initial chromosomes of a generation, (ii) Apply a small number of crossovers, (iii) Use a selection procedure similar 
to that of Algorithm [5J Nevertheless, as we have discussed, it is not known how to construct the procedure (ii) as 
a unitary operation without the expensive process of looking up classical data of chromosomes in the generation. 
Therefore, it is the fact that a meaningful quantum counterpart is not easily found for a common classical genetic 
algorithm. 

In view of the search space covered by each generation, a simultaneous crossover is, of course, desirable. Using 
a superposition for this purpose was discussed [Sj, but not developed previously. In this sense, we have made a 
meaningful improvement by introducing Algorithmic 

There have not been many studies on quantum genetic algorithms in the literal meaning (i.e., not in the meaning 
of quantum-inspired classical algorithms) so far. It is hoped that useful quantum extensions other than our proposal 
will be developed in future. 



6 Summary 

We have proposed a semiclassical quantum genetic algorithm whose crossover, mutation, and selection procedures 
have been constructed as quantum routines so that quantum parallelism is effectively used. Our complexity analysis 
has shown that the algorithm achieves a quadratic speedup over its classical counterpart. 
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Appendix: An example of constructing the pseudo-randomizer R 

In this appendix, we show an example to construct the pseudo-randomizer R used in Algorithms [2l [H and|4l It should 
map a c-bit string to a pseudo-random n-bit string except for Oq • • • Oc-i that is mapped to Oq • • • 0„_i. Its internal 
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circuit complexity should be poly(cn). As we use a quantum circuit to realize it in Algorithm [21 it is desirable to 
employ a circuit structure that is originally unitary. 

Consider inputs a e {0, 1}"^. We design a circuit that maps ao • • • ac-i% ■ ■ ■ 0„_i to ao • ■ ■ ac-iKo • ■ ■ Kn-i with 
K — R{a) an n-bit pseudo-random number (here, oq ■ • ■ flc-i and kq ■ • • ^n-i are the binary representations of a and k). 
By the definition of i?, the circuit preserves Oq ■ • • Oc-iOq ■ • • 0„_i. This circuit is generated by the following function. 

FUNCTION gen_r_circ(): 

Comment: We use wires vq, . . . , Vc-i^wq, . . . , w„_i. 
FOR ii — TO c- 1: 

(1) Use a random number generator to generate an 7i-bit integer 7. Write its binary representation as 
70 • • -Tn-i- 

(2) Using the wire Vi as the control wire (namely, the control bit), output the gate controlled-{XQ)^" (E) 
• • • (X) {Xn-i)^"~'^ with Xk the bit flip gate acting on the wire Wk {k — 0, . . . ,n — 1). In this gate, the 
bit flips are active under the condition that Vi = 1. 

END FOR 

Obviously, the circuit output from this function can be directly used as a quantum circuit. Using the circuit C = 

gen_r_circ(), we have |a)|0) \a)\R{a)) Va e {0,1}'^. The circuit complexity of C is 0{cn) because, for each 
I, at most n CNOT gates are used to decompose the gate output from step (2). In addition, geii_r_circ() spends 
0{c poly(?i)) time when a common random number generator [Hlj l is used in step (1). 



Note that the function gen_r_circ() is called only once for each t, in the beginning of step 1 in Algorithms 
[21 [31 and[H We have only to reuse the circuit C for the use of the pseudo-randomizer until t is incremented. 

It is expected that outputs from the circuit C possess good uniformity if we use a good random number generator 
in step (1) of gen_r_circ() for generating C. Let us write 7 as 7(1) to emphasize its dependence on i. For a nonzero 
input uq ■ ■ ■ Oc-i, the fcth bit of the output R(a) is X]i=o '^i ' 7fe(*) mod 2. This indicates that, for two different inputs 
a and a', the kth bits of R{a) and R{a') differ with probability 1/2 in the ideal case where 7(«)'s are generated from 
a true random number generator. This is because a and a' differ by at least a single bit. It also indicates that two 
different bits, the kth and the fc'th bits, of R{a) for a nonzero input a differ with the probability 1/2 in the ideal case. 
This is because jk{i) and 7fc'(i) differ with the probability 1/2. 

Now we show the result of our numerical test of C. We tried statistical tests of randomness [Hol [s^ to test pseudo- 
random numbers output from C, using NIST's Statistical Test Suite (STS) (version 2.1.1) [521. We set c = 10 and 
n = 32. Mersenne Twister (MT) (version mtl9937ar) 51] was used to generate 7 in step (1) of gen_r_circ(). We used 
the seed value 121212 and did not reset MT during the circuit generation. The circuit C output from gen_r_circ(), 
of course, consisted of 10 outputs from step (2). For this C, we used the inputs a £ {0, l}'^\{0o • • • Oc-i} from smaller 
to larger and obtained corresponding outputs R{a) by numerical computation. We obtained 1023 x 32 bits in total 
in the outputs, since 2^" — 1 = 1023. We regarded them as a serial bit string from left to right and used STS in its 
default setting to test the string. In the execution of STS, we used 25 binary sequences with length 1200 as samples 
from the string. The following tests were tried with the default parameter values in STS: Frequency Test, Block 
Frequency Test, Cumulative Sums Test, Runs Test, Longest-Run-of-Ones Test, Binary Matrix Rank Test, Spectral 
DFT Test, and Serial Test. The string passed the tests except for Binary Matrix Rank Test. It should be noted 
that the input length was too small for the binary matrix rank test [53| . In addition, randomness is not very strictly 
required for the use in evolutionary computing. Therefore, considering the tests that the string passed, we may claim 
that gen_r_circ() generates a usable pseudo-randomizer circuit for our algorithm. 

We conducted another test: We generated ten circuits by calling gen_r_circ() ten times without resetting MT, 
using the seed value 676767. For each circuit, we performed the same process as above to obtain the serial bit string. 
We obtained ten serial bit strings in total and tested the concatenated string using STS. As samples input to STS, 
we used 25 binary sequences with length 12000. The concatenated string passed the tests except for Binary Matrix 
Rank Test and Spectral DFT Test. It was unexpected that it did not pass the spectral DFT test. It requires a further 
investigation to reveal the reason of this phenomenon. 

The results of the first and the second tests are summarized in Table [U In summary for this appendix, we 
found that a pseudo-randomizer circuit whose outputs possess enough randomness can be generated by the function 
gen_r_circ() for the use in our algorithm. It is hoped that the function will be improved so as to achieve better 
randomness for the sake of general use. 
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